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Abstract 

Acoustic wave generation by turbulence in the stratified, moist atmosphere is studied. It is 
shown that in the saturated moist air turbulence in addition to the Lighthill's quadrupole and the 
dipole sources of sound related to stratification and temperature fluctuations, there exist monopole 
sources related to heat and mass production during the condensation of moisture. We determine 
acoustic power of these monopole sources. Performed analysis shows that the monopole radiation 
is dominant for typical parameters of strong convective storms. Obtained results are in good 
qualitative agreement with the main observed characteristics of infrasound radiation by strong 
convective storms such as total acoustic power and characteristic frequency. 

PACS numbers: 47.27.Sd,43.28.+h 
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I. INTRODUCTION 



It is long known that strong convective storms, such as supercell thunderstorms are 
powerful sources of infrasound IlL M. Detailed observations of convective storm generated 
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infrasound [1|, |3j provides that at least two different group of infrasonic signals could be 
identified. The first group, with characteristic period about 1 s, has been found to be in 
strong connection with proto-tornadic structures, funnel clouds and tornadoes. Based on 
coincident radar measurements of tornados, which show strong relationship between funnel 
diameter and infrasound frequency, it is usually supposed that these infrasound waves are 
generated by radial vibrations of the vortex core 4j . The second group of infrasound signals 
has the periods from 2 to 60 s. Usually the emission appears about 1 hour before observation 
of tornado. Although detection of these waves are known to be strongly correlated with 
formation of tornado |5j| , it is usually supposed that these waves are not related with tornado 
itself and are caused by convective processes that precede tornado formation. The acoustic 
power radiated by convective storm system could be as high as 10 7 watts 5|. Although 
several reasonable mechanisms have been suggested to explain this acoustic radiation, the 
physical mechanism of the process remains unexplained % |5|, |6j. 

Broad and smooth spectrum of the observed infrasound radiation indicates that tur- 
bulence is one of the promising sources of the radiation. Lighthill's acoustic analogy 
represents the basis for understanding of the sound generation by turbulent flows. Valida- 
tion of this theory has been shown by various experiments and numerical simulations (see, 



e.g., Refs 



. 0,y. 
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llj and references therein). In this approach the flow is assumed to 



be known and the sound field is calculated as a small by-product of the flow. According 
to this theory in the case of uniform background thermodynamic parameters interaction of 
turbulent vortices provides quadrupole source of sound. The acoustic power of the source 



was estimated by Proudman 12]. But usage of this estimation for the infrasound radiation 
from convective storms usually leads to the underestimation of the acoustic power 
It requires characteristic velocity of the turbulence to be much greater than it exists in 
any terrestrial storm system. Recent analysis of a non-supercell tornado storm simulation 

n 

performed by Nicholls, Pielke and Bedard [14( suggests that the occurrence of the high fre- 
quency infrasound coincided with the development of considerable small-scale turbulence 
that may have caused small-scale latent heating fluctuations which appear to be the main 



2 



mechanism responsible for generating the infrasound in this simulation. 

In the presented paper we study acoustic radiation from turbulent convection taking into 
account effects of stratification, temperature fluctuations and moisture of the air and using 
Lighthill's acoustic analogy. Formulation of the generalized acoustic analogy [15] implies: 
(i) dividing the flow variables into their mean and fluctuating parts; (ii) subtracting out the 
equation for the mean flow; (iii) collecting all the linear terms on one side of equations and the 
nonlinear terms on the other side; (iv) treating the latter terms as the known terms of sound. 
We show that in the saturated moist air turbulence in addition to the Lighthill's quadrupole 
and known dipole sources of sound related to stratification and temperature fluctuations, 
there exist monopole sources related to heat and mass production during the condensation of 
moisture. It appears that infrasound radiation from convective storms should be dominated 
by acoustic source related to the monopole sources related to the moisture of the air. We 
show that for the typical parameters of the strong convective storm the acoustic output of 
this monopole source is two orders of magnitude stronger than Lighthill's quadrupole source, 
whereas the dipole radiation related to temperature inhomogeneities is of the same order 
as radiation of Lighthill's quadrupole source. The dipole source related to stratification 
and the dipole and quadrupole sources related to inhomogeneity of background velocity are 
inefficient sources of sound. The total power of the source related to moisture is of order 
10 7 watts for the typical parameters of the strong convective storms, in qualitative agreement 
with observations j^, 3, 6]. 

The paper is organized as follows: In Sec. [Til equations governing sound generation 
by turbulence for moist atmosphere are obtained in the framework of Lighthill's acoustic 
analogy. Various sources of acoustic radiation are analyzed in Sec. IIHI Application of the 
obtained results to infrasound generation in strong convective storms is made in Sec. IIVI 
Conclusions are given in Sec. [V] 

II. GENERAL FORMALISM 

The dynamics of convective motion of moist air is governed by continuity, Euler, heat, 
humidity and ideal gas state equations: 

^ + PV.V = 0, (1) 
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p— + 2pft x v = -Vp - pV$, (2) 

J—J L 

1 Dt~ V DV {d} 

= P 1 = P 1 (A) 

P RTl-q + q/e RTl + aq U 

where v, 2f2xv, p andp are velocity, Coriolis acceleration, density and pressure respectively; 
D/Dt = d/dt + v • V is Lagrangian time derivative; L v is the latent heat of condensation 
and q is the mass mixing ratio of water vapor (humidity mixing ratio) 

q = 7' (5) 

where p u is the mass of water vapor in the unit volume; e = m v /md ~ 0.622 is the ratio of 
molecular masses of water and air; a = 0.608 and R is the universal gas constant. 

In the set of Eqs. ©-(jlj) diffusion and viscosity effects are neglected due to the fact 
that they have minor influence on low frequency acoustic wave generation as well as its 
propagation. 

In the future analysis we also assume Q = 0. As it is well known 16], Coriolis effects 
are negligible for mesoscale convective system dynamics. On the other hand, when the 
frequency of acoustic waves Q a satisfy the condition Q a ^> Q, Coriolis effects have also 
negligible influence on acoustic wave dynamics. 

The main idea of Lighthill's acoustic analogy is reformulation of the governing equations 
in the form suitable for the study of acoustic wave radiation process. To proceed in this 
direction one have to choose appropriate "acoustic variable", that describes acoustic waves 
in the irrotational regions Generalized Bernoulli's theorem fl suggests that 

the total enthalpy 

p v 2 

B = E+ - + — + $, (6) 
P ^ 

where $ is gravitational potential energy per unit mass, E is internal energy and V$ = — g, 



is one of the possible appropriate choices 18]. B is constant in the steady irrotational flow 



and at large distances from acoustic sources perturbations of B represent acoustic waves. 

For derivation of acoustic analogy equation in terms of the total enthalpy it is useful to 
rewrite Euler's equation in the Crocco's form 

p^ + VB = -uxv + TVs, (7) 
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where u is vorticity, T is temperature, s is specific entropy and 

Tds = dE + pd(^ =dB-^-d$-d(^J . (8) 
From the thermodynamic identity 

*P=( d A W^l W 



v^P/.,, \ 9s /p,g \ d 1ys, P 

where the subscripts serve as the reminders of the variables held constant, using Eqs. PJ 
we obtain 

dp = —dp — —ds -\ — dq, (10) 

cj c p 1 + aq 



where 



dp 



s,q 



is the sound velocity and 



^T[%) (12) 



<9T 



is the specific heat of the air. 

Eliminating convective derivative of the density from Eq. (pQ) using Eq. (flUj) we have 

1 Dp ^ 1 Ds 1 Dq 

-^7— + V-v = + -. (13) 

pc^ L>t c p Dt 1 + aqDt K J 



Subtracting the divergence of Eq. flxj) from time derivative of Eq. (flBl and using Eq (Ej) 
after long but straightforward calculations we obtain 

D (1 D\ V P -V 2 \ 

S-x + 5 r + 5, + 5^ + 5 7 , (14) 
where 7 = c p / c v is the ratio of specific heats and 

5l "( v+ S) >xv) ' (l5) 

S T = - (V + • (TV,) , (16) 
9 f 1 TDs\ , ^ (T Ds\ 



s "^Ui«J +(v ' v) UW' (17) 
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^-^(^(W)p-^(vV)g). (19) 



#7 (dq dp 

Eq. ffT4j) is suitable for identification of different acoustic sources and the study of their 
acoustic output. 

The nonlinear wave operator on the left of Eq. ( fl4l) is identical with that governing the 
propagation of sound in irrotational, homentropic flow. Therefore the terms on the right 
may be identified as acoustic sources. Propagation of infrasound in the atmosphere was 
intensively studied by different authors (see, e.g., Ref. [3] and references therein) and will 
not be considered in this paper. 

For simplification of further analysis of the acoustic output of different sources we make 
several standard assumptions: 

(a) Studding acoustic wave generation process for low Mach number flow, all the connec- 
tive derivatives in Eq. ffl4|) can be replaced by time derivatives d/dt j^; 

(b) For acoustic waves with the wavelength A not exceeding the stratification length scale 



tit 



\<H = -^^ 10 4 m, (20) 
one can also neglect the influence of stratification on the acoustic wave generation process 



and consider background thermodynamic parameters in Eq. (I14p as constants 20]. 

(c) Neglecting nonlinear effects of acoustic wave propagation and scattering of sound by 
vorticity and taking into account that M <§C 1, for the acoustic pressure in the far field we 
have 

p , (x,t)^po5(x,t). (21) 

(d) Eq. ffT4l) is equivalent to initial set of Eqs. CQ)-(JI]) and therefore it describes not only 
acoustic waves, but also the instability wave solutions that are usually associated with large 
scale turbulent structures and continuous spectrum solutions related to "fine-grained" tur- 



bulent motions 
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2l| . In the presence of any kind of inhomogeneity, such as stratification 
or velocity shear, linear coupling between these perturbations is possible, and in principle 
acoustic waves can be generated by both instability waves and continuous spectrum pertur- 
bations. But in the case of low Mach number (M -C 1) flows both kind of perturbations are 
very inefficient sources of sound. The acoustic power is proportional to e~ l l 2M2 and e~ n5 ^ 2M 



for instability waves and continuous spectrum perturbations respectively 22j. In the last 
expression 5 is the ratio of length scales of energy containing vortices and background ve- 
locity inhomogeneity (V/d z V). In the case of supercell thunderstorm M ~ 0.1 — 0.15 and 



5 ~ 10~ 2 , therefore both linear mechanisms have negligible acoustic output and attention 
should be payed to sources of sound related to nonlinear terms and entropy fluctuations that 
will be studied below. 

With these assumptions Eq. (|14p simplifies and reduces to 

I fl d 2 



with 



jH) y^^-^ 2 )p' = S L + S T + S, + S q + S m , (22) 



S L RiV-(wxv), (23) 
S T ~ -V ■ (TVs) , (24) 

^-'Sg^'-IH- (25) 



d 2 q 
1 + aqdt 2 



First three terms on the right hand side of Eq. (1221) represent well known sources of 
sound: the first term represents Lighthill's quadrupole source [7]; the second term is dipole 
source related to temperature fluctuations [23j; S 7 is monopole source related to variability of 
adiabatic index, that usually have negligible acoustic output [18J and will not be considered 
in the presented paper; Eq. (122]) shows that in the case of saturated moist air turbulence 
there exist two additional sources of sound. S q and S m are monopole sources related to 
nonstationary heat and mass production during the condensation of moisture, respectively. 



III. ANALYSIS OF DIFFERENT SOURCES 



For estimation of different acoustic sources we follow the standard lg, |23j procedure. 
Namely, using the wave equation free space Green function 

5{t-t' - |x-x'|/c s ) 



G < M ' X ' X,) 4^|x-x'| • 

acoustic pressure fluctuations corresponding to a source Si can be written as 



(28) 



Pi{x,t) 



1 f [St]t=t * dV, 



4:71 C 2 J |x — X' 



(29) 



where t* = t — |x — x'l /c. 



Calculating acoustic radiation in the far field (|x| ^> |x'|), we can use following expansions 

(30) 



x • x 



s, (t> )„ Mt _M) + ^»W t _M,, pi) 



and using plane wave approximation for far field derivatives 

d Xi d 

r^j 

dxi c a \x\dt' > 

for the Lighthill's source we obtain 



(32) 



M'A = -{^W (33) 

which corresponds to quadrupolar radiation field. 

p' L can be estimated in terms of characteristic velocity v and length scale I of energy 
containing turbulent eddies. Fluctuations in ViVj in different regions of the turbulent flow 
separated by distances greater then / tend to be statistically independent, and therefore 
generation of sound can be considered as by a collection of F/l 3 independent eddies, where 
F is the volume occupied by the turbulence. The dominant frequency of the motion ~ v/l, 
so the wavelength of the radiated sound Q ~ l/M t , where M t = v/c s <C 1 is turbulent Mach 
number. Therefore, each eddy is acoustically compact. Acoustic pressure generated by 
single eddy is p' L1 ~ (l/\x\)pov 2 M 2 , and acoustic power it radiates Nn ~ 4n |x| 2 pf 1 / poc s ~ 
Pov 3 l 2 M 5 , that corresponds to Lighthill's eighth power law. For total acoustic power this 
yields Praudman's estimate 12] 

N L ~ y-y-P 1 - (34) 
lei 

Similar arguments can be used for estimation of acoustic power of thermo-acoustical 
source Sj_ related to density (and therefore temperature) fluctuations, that produce dipole 
source 18(. The physics of this kind of acoustic radiation is the following: "hot spots" or "en- 
tropy inhomogeneities" behave as scattering centers at which dynamic pressure fluctuations 
are converted directly into sound. The acoustic power is 

pATV AT 2 . . 

met M 2 T 2 1 ' 

where AT denotes the rms of temperature fluctuations. 
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Acoustic sources S q and S m are related to the moisture of the air. They produce monopole 
radiation and physically have the following nature: suppose there exist two saturated air 
parcels of unit mass with different temperatures 7\ and T 2 and water masses m u (Ti) and 
m v (T 2 ). Mixing of these parcels leads to the condensation of water due to the fact that 

2m„(T 1 /2 + T 2 /2) < m^Tj + m u (T 2 ). (36) 

Condensation of the water leads to two effects, important for sound generation: pro- 
duction of heat and decrease of the gas mass. Both of these effects are known to produce 



monopole radiation [18l. l23l|. Consequently, turbulent mixing of saturated air with different 
temperatures will lead not only to the dipole thermo-acoustical radiation fl35|) . but also to 
the monopole radiation. 

According to Eqs. ([5]) and ( 1361) for humidity mixing ratio fluctuation q' s we have 

q' s = q s (T + T') + q s (T - T') - 2q s (T) . (37) 

In the limit T'/T <C 1 this yields 

fli - (38) 
Substituting (JXTJ) into p| and using ([3"8]). pOl) and (13T1) we obtain 

PiCM) = -^Sfe^l^ / T'(x',t)T>(x>,t)d^, (39) 



4ttc2|x| dT 2 dt 2 _ 

which corresponds to monopole radiation field. 

For total acoustic power radiated by monopole source related to the moisture we have 

47r|x 



iV g = ^L(p'(x > t) J /(x > t))~ (40) 
Poc s 



P0l2Ll (§fi) % j d 3 x'd 3 x"(T\x',t)TXx>,t)T>(x",t)T'(x",t)) 



Fluctuations of temperature in different regions of the turbulent flow separated by distances 
greater then length scale I of energy containing eddies are not correlated and therefore 
the integral in Eq. (T4T|) can be estimated as Fi/ 3 AT 4 , where F\ is the volume occupied by 
saturated moist air turbulence. Taking also into account that the characteristic timescale of 
the process is turn over time of energy containing turbulent eddies l/v finally obtain 

Pol 2 LlAq 2 M? „ _ 7 2 LjA_g 2 F x 

ft4 



q lc„ Mfct F K J 



where Aq is the rms of humidity mixing ratio perturbations, the acoustic power of the 
source related to the gas mass production we obtain 



p a 2 c 3 s Aq 2 M? a 2 Aq 2 F 1 M 



I 



M 4 F 



(42) 



IV. APPLICATION TO INFRASOUND GENERATION BY STRONG CONVEC- 
TIVE STORMS 

In this section we apply our findings to study infrasound generation by strong convective 
storms. Taking for typical parameters of supercell storms v ~ 5 m/s, AT ~ 3° K ,T = 

270° K and c s = 330 m/s and using Eqs. (151]) and (135]) we see that the dipole radiation related 
to temperature inhomogeneities is of the same order as radiation of LighthiU's quadrupole 
source. 

Combining Eqs. (IHl-(jl2l). using L u ps 2.5 x 10 6 m 2 /s 2 and 7 « 1.4 we obtain 



N m 



10 3 



(43) 



therefore acoustic power of the source related to the gas mass production is negligible com- 
pared to the radiation related to the heat production. 

Estimation of Aq is a bit more difficult. For saturation specific humidity we use Bolton's 



formula 



24) 



3.8 

q s « — exp 
Po 



17.67T C 
T c + 243.5 



(44) 



where T c = T — 273.15 is the temperature in degree Celsius and po is atmospheric pressure 
in mb. Taking into account (1381 and using p w 800 mb, we obtain 



Aq 



6.8 • 10 4 



exp 



17.67T, 



AT 2 = f(T c 



AT 2 

J^2 



(45) 



(243.5 + T c ) 4 ~" r \T C + 243.5 / 

Note, that due to the numerator in the exponent f(T c ) strongly depends on temperature, 
e.g., /(T c = 10°)//(T o = 0°) « 2. 
Using Eqs. (151]) and (jlT]) we obtain 

2 / 2 (T c ). (46) 







2 


" AT " 


4 


"iV 





c 2 








T 



For our analysis we assume T c = 0°C, (/(0) ~ 1.66) and T ~ 125 km 3 [6]. For estimation 
of iq we note that that for atmospheric convection saturation level is usually at the height 



10 



m 1 — 1.5 km. Taking also into account that f(T c ) rapidly drops with the decrease of T c , 
one can expect that main acoustic radiation will be produced at the heights (1.5 — 4) km, 
and consequently we assume Fx ~ 0.5F. Then Eq. ( I4"6"i) yields 

|- 2 x 10 2 . (47) 

Therefore, we conclude that infrasound radiation of supercell storm should be dominated by 
the monopole source related to the heat production during water condensation. Assuming 
additionally constant of proportionality in Eq. ( 1341) equal to 100 [f], [3| and I « 10 m for 
total power of the radiation we obtain 

N q w 2.4 x 10 7 watts, (48) 

in qualitative agreement with observations 

As it was mentioned above, the characteristic frequency of the emitted acoustic waves 
Q ~ v/l, and using the characteristic values of the velocity and length scale we obtain for 
the period r ~ 10 s. 



V. CONCLUSIONS 



In the presented paper we have considered acoustic radiation from turbulent convection 
in the framework of generalized acoustic analogy taking into account effects of stratification, 
temperature fluctuations and moisture of the air. Analysis shows existence of monopole 
sources related to heat and mass production during the condensation of moisture in the 
saturated moist air turbulence, in addition to the Lighthill's quadrupole and known dipole 
sources of sound related to stratification and temperature fluctuations. It has been shown 
that for the typical parameters of the strong convective storms infrasound radiation should 
be dominated by monopole source related to the moisture of the air. The total power of the 
source related to moisture is of order 10 7 watts, in qualitative agreement with observations 
of strong convective storms 
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